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We demonstrate how to separate the rotational degrees of freedom in 
a quantum iV-body problem completely from the internal ones. It is shown 
that any common eigenfunction of the total orbital angular momentum (£) 
and the parity in the system can be expanded with respect to (21 + 1) 
base- functions, where the coefficients are the functions of the internal vari- 
ables. We establish explicitly the equations for those functions, called the 
generalized radial equations, which are (2£ + 1) coupled partial differential 
equations containing only (3N — 6) internal variables. 

PACS number(s): ll.30.-j, 03.65.Ge, and 03.65.Fd 

Symmetry is an important property of a physical system. The symmetry of a quantum 
system can simplify its Schrodinger equation and remove some variables in the equation. 
The simplest example is the hydrogen atom problem, where, due to the spherical sym- 
metry, the wavefunction is expressed as a product of a radial function and a spherical 
harmonic function, 



and the Schrodinger equation reduces to a radial equation with only one radial variable. 
For a quantum iV-body problem with a pair potential, the Schrodinger equation is in- 
variant under the spatial translation, rotation, and inversion. From those symmetries, 
one should be able to separate the motion of center-of-mass and the global rotation of 
the system from the internal motions so as to reduce the Schrodinger equation to the 
generalized "radial" equation that contains only internal variables. However, this prob- 
lem has not been solved. In this letter we will solve this problem completely. Using the 
appropriately chosen (3iV — 6) internal variables and the (2£ + 1) base-functions for the 
total orbital angular momentum £, we establish explicitly the generalized radial equations 
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without any approximation. Only (3iV — 6) internal variables are involved in both the 
generalized radial functions and the equations. 

Denote by rj the position vectors of iV-particles with masses m 3 - in the laboratory 
frame (LF), respectively. The Schrodinger equation for the iV-body problem is 

N 

— (h 2 /2) m f A r, * + = EV, (2) 

3=1 

where V is assumed to be a pair potential, depending on the distances of each pair of 
particles. Therefore, the potential V is a function of only the internal variables. It is well 
known that, due to the translation symmetry of the system, the motion of center-of-mass 
can be separated completely from others by making use of the Jacobi coordinate vectors 
in the center-of-mass frame (CF) [1-3], 

where Wj = Y^t=j m t- I n CF, the Laplace operator and the total orbital angular momen- 
tum operator L can be directly expressed with respect to Rfc: 

N N-1 

A = E <A rj = £ A Rfc , 

3=1 k=l 

N N-1 

L = -ihJ2 r i x Vr, = -ih R fe x VR t , (4) 

3=1 k=l 

The Laplace operator obviously has the symmetry of 0(3N — 3) group with respect to 
(3iV — 3) variables. The 0(3iV — 3) group contains a subgroup SO (3) x 0(N — 1), where 
SO (3) is the usual rotational group. The space inversion and the different definitions 
for the Jacobi vectors in the so-called "Jacobi tree" [4] can be obtained by 0(N — 1) 
transformations. For the system of identical particles, the permutation group among 
particles is also a subgroup of 0(N — 1) group. 

Because of the spherical symmetry, the angular momentum is conserved. The hydrogen 
atom problem is a typical quantum two-body problem, where there is only one Jacobi 
coordinate vector, usually called the relative position vector r. For a quantum iV-body 
problem, equation (1) should be generalized in three aspects. The first is how to define 
the internal variables, which describe the internal motions completely. The second is 
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how to find the complete set of the independent base-functions with the given angular 
momentum. The total wavefunction is expanded with respect to the base-functions, where 
the coefficients are the generalized radial functions which only depend on the internal 
variables. The last is how to derive the generalized radial equations that only contain 
(3iV — 6) internal variables. As a matter of fact, these three aspects are connected. The 
parity should also be considered in the generalization. Due to the spherical symmetry, one 
only needs to study the eigenfunctions of angular momentum with the largest eigenvalue 
of L3 (m = £), which are simply called the wavefunctions with the angular momentum 
£ in this letter for simplicity. Their partners with the smaller eigenvalues of L3 can be 
calculated from them by the lowering operator L_. 

Denote by R = R(a, (3, 7) a spatial rotation, transforming CF to the body-fixed frame 
(BF), and by £ all the internal variables in a quantum iV-body problem for simplicity. 
Although Wigner did not separate the motion of center-of-mass by the Jacobi vectors, he 
proved from the group theory that any wavefunction with the angular momentum £ in 
the system can be expressed as follows (see Eq. (19.6) in [5]): 

*fta,/3,7,0= E /3,7)W0, (5) 

q=-t 

where we adopt the commonly used form of the D-function [6]. In Eq. (5) Dg q (a, (3, 7)* 
plays the role of the base-function with the angular momentum £, and ip q (£) is the gener- 
alized radial function. What Wigner proved is that there are only (2£ + 1) independent 
base-functions with the angular momentum £. Unfortunately, due to the singularity of 
the Euler angles, the generalized radial equations are very difficult to derive based on 
Eq. (5). Wigner did not discuss the generalized radial equations, and to our knowledge, 
those equations have not yet been established in the literature. It is obvious that the 
generalized radial equations are very easy to obtain for the S wave [7]. However, it seems 
quite difficult to obtain even for P wave in a three-body problem [8,9]. 

Recently, a coupled angular momentum basis was used to prediagonalize the kinetic 
energy operator [10], where some off-diagonal elements remain non- vanishing. In their 
calculation, the function with a given angular momentum was combined from the partial 
angular momentum states by the Clebsch-Gordan coefficients. Since the partial angular 
momenta are not conserved, one has to deal with, in principle, an infinite number of 
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the partial angular momentum states. This problem also occurs in the hyperspherical 
harmonic function method and its improved ones [2,4,11,12]. However, as Wigner proved, 
only (2£ + 1) partial angular momentum states are involved in constructing the base- 
functions with the angular momentum £. 

Arbitrarily choose two Jacobi coordinate vectors, say Ri and R 2 . Let Ri be parallel 
with the Z-axis of BF, and R 2 be located in the XZ plane with a non-negative X- 
component in BF. The rotational degrees of freedom of the system are described by a 
rotation R(a,(3, 7), transforming CF to BF. Define (3iV — 6) internal variables, which 
should be invariant in the global rotation R(a,j3, 7): 

= Rj-Ri, ^ = Rj-R 2 , = Rj-(RiAR 2 ), 1<J<(N-1), (6) 

where 7/1 = £ 2 and Ci = C2 = 0. It is worth mentioning that ^ and r/j have even parity, 
but Q has odd parity. From them we have 

Slj = (Ri A Rj) ■ (Ri A R 2 ) = faj - 

Uj = (R 2 A Rj) ■ (Ri A R 2 ) = &Vj - ViZj, 

Rj ■ Rk = &2 1 (tyVk — ^jik + (j(k) , (7) 

where Vti = uj 2 = 0, and Vt 2 = —uj\ = (Ri A R 2 ) 2 . 

Recall that two Jacobi vectors Ri and R 2 completely determine BF and three Euler 
angles. The base-functions with the angular momentum £ should be combined from the 
products of two spherical harmonic functions Y^(Ri) and Y^,(R 2 ) by the Clebsch-Gordan 
coefficients (q,m,p,m'\£,(m + m')). Define [3,11] 

Qf (Ri, R 2 ) = (^ + ^i 2 )^(^i + ^ 22 )^ {{Ru + iRu)R23 _ Ri3{R2i + iRm)y 

r<q<£, r = 0,l. (8) 

where Rj a is the ath component of the Jacobi vector R^-. Q £ q T (R 1 ,R 2 ) is the common 
eigenfunction of L 2 , L 3 , A Rfe , and the parity with the eigenvalues £(£ + 1), £, 0, and 
( — l) i+T , respectively. As a matter of fact, the following combination of products of two 
spherical harmonic functions is proportional to Q e q T (Ri, R 2 ) 

e ^v^R^r^^ (9) 

m 
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where C is a normalization factor. Now, we come to the theorem. 



Theorem. Any function ^ A (R l5 • • • , Rjy_i) with the angular momentum £ and the parity 
(— 1) £+A in a quantum iV-body problem can be expanded with respect to Q e q T (Ri, R 2 ) with 
the coefficients ^^(£,1], (), which depend on (3N — 6) internal variables: 



*?(Ri, • • • , Rjv-i) = EE V, C)QT( R i» R 2)' 

T =o g=-r 

V'Jr (f » *7> = i>qr (6, ' ' ' , ^2, • • • , VN-1, Cs, ' " " , Ov-l), 



(10) 



where the parity of t/>^(£, 77, C) is ( — 1 



a-r 



Equation (5) coincides with Eq. (10), because either of the set of D\(a, (3, 7)* and the 
set of Q e q T (Ri, R 2 ) is a complete set of base-functions of the angular momentum. However, 
equation (10) has three important characteristics, which make it possible to derive the 
generalized radial equations. The first is that the Euler angles do not appear explicitly 
in the base-functions Qg T (Ri,R 2 ). The second is the well chosen internal variables (6). 
The third is that the internal variables Q have odd parity. It is due to the existence 
of Q that the base-functions (3g°(Ri, R2) and Qg 1 (Ri,R 2 ) appear together in one total 
wavefunction. By comparison, all the internal variables in a quantum three-body problem 
have even parity (Q = 0) so that in a total wavefunction with a given parity, only the 
base-functions with the same parity appear [5,12]. 

Now, substituting Eq. (10) into the Schrodinger equation (2) with the Laplace operator 
(4), we obtain the generalized radial equations by a straightforward calculation: 



Atf* + 4 {g% + (£- q)d m }^ + 2qd^_ 1)0 + 2(£ - g)%< +1)0 



1, ex 



ex 



i.ex 



N-l 



1 ex 

J q0 



+ £ { [-Ujgd^ + Qjit - q)d Vj + rj 2 Qqd Cj + ^Q(£ - q)d Cj ] ifi 

i=3 

-q [iOjd Vi + bQdcs] ^ A _ 1)0 + (t~q) [toAj - tf}+i)o 
-irj 2 q(q - 1) [Qd Vj - Q^.] 
-iq [mQqdt, - U 3 {2£ -2q+ l)d Vj + m^qd Cj + ^{2£ - 2q + l)d Cj ] ^ 
- q) [60(2? + 1)% - ^Q(£ - q)d Vj + 6^-(2g + l)d Cj + ^(£ - q)d Cj ] 

-qW-q-l) [CA> + "A] 1^2)1} = " {2/h 2 ) [E - V] ^g, (11a) 



N% + 4 {g% + (£ - g + 1 A> ^ + 2(g - + 2{l - 



i, ex 



N-l 



+ £ 2Q 2 " 1 { [-c^. + n,-(€ - 9 + 1 A + 7720^ c . + " q + 1 A] 



i=3 



CA - n A ^ 



/A 



j I ^(9-1)0 



CA 



W A 



^g} = -(2/n 2 )[s-v]^g, (nfc) 

A^(C, C) = {46^ + 4^ + (6 + 772) + 46 (% + % + 6 («% + 9, 2 ) 

JV-1 

+ E Mi + ^ + (vM - fa + c|) 0J + aj,) 

•5V (Hi + + "J + 6C| + V2C-) dl + 4 (i/'A + CA) <% 

+4 + CA) d m + 2 (vA, + CA) % + 26^A] } ^ 0- (He) 

Due to the limited size of a letter, we have to leave the proof of the theorem and the 
detailed calculation elsewhere. When establishing BF we arbitrarily choose two Jacobi 
coordinate vectors Ri and R2. Those two vectors may be replaced with any other two 
Jacobi vectors. One may change the choice according to the characteristics of the physical 
problem under study, such as some or all particles in the quantum iV-body problem are 
the identical particles. 

In deriving the generalized radial equations, the key is to discover the base-functions 
Qg T (Ri,R2) of the angular momentum and to choose the right internal variables, some 
of which have odd parity. From Eq. (9) we see that only finite number of partial angular 
momentum states are involved in constructing the base- functions Qg T (Ri,R2)- Namely, 
the contributions from the remaining partial angular momentum states have been incor- 
porated into those from the generalized radial functions. 

The two features in this method, that the numbers of both functions ip£(€,r},£) and 
equations are finite, and they depend only on (3iV — 6) internal variables, are important for 
calculating the energy levels and wavefunctions in a quantum iV-body problem. In fact, 
in the numerical experiment for the quantum three-body problem by the series expansion, 
much fewer terms have to be taken to achieve the same precision of energy than with other 
methods. The calculation error will be less in comparison with the method to truncate 
the series on the partial angular momentum states. As the number of the particles in the 
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system increases, we believe, to remove three independent variables will greatly decrease 
the calculation capacity requirement. 
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